Formation of multi-quanta vortices in superconductors: 
electronic, calorimetric and magnetic evidence 



Ben Xu, 1 M. V. Milosevic, 1 Shi-Hsin Lin, 2 ' 1 F. M. Peeters, 1 and B. Janko 2 

1 Departement Fysica, Universiteit Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium 
^Department of Physics, University of Notre Dame, Notre Dame, IN 46556, USA 

(Dated: July 4, 2011) 

The ground state with vorticity larger than one in mesoscopic superconductors in applied magnetic 
field may manifest as a 'giant'-vortex, where all vortices coalesce into a single singularity of the 
order parameter. Such a multi-quanta vortex may split into individual vortices (and vice versa) as a 
function of e.g. applied current, magnetic field or temperature. Here we show that such transitions 
can be identified by heat-capacity measurements, as the formation or splitting of a giant-vortex 
results in a clear jump in measured heat capacity vs. external drive. We attribute this phenomenon 
to an abrupt change in the density of states of the quasiparticle excitations in the vortex core(s), 
and further link it to a sharp change of the magnetic susceptibility at the transition - proving that 
formation of a giant-vortex can also be detected by conventional magnetometry. 
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The influence of quantum confinement on supercon- 
ducting condensates is certainly one of the prominent 
research directions in low-temperature physics of the last 
decade. Besides the ever intriguing properties of high- 
temperature superconductivity, vortex matter in conven- 
tional but mesoscopically tailored superconducting sam- 
ples has generated tremendous interest and activity in 
the wide scientific community. Indeed, vortex configu- 
rations in confined condensates are of direct relevance 
to cold gases and Bose- Einstein condensates [l|; interac- 
tion of vortices with artificial pinning sites has analogues 
in various colloidal systems and molecular crystals 
the inhomogeneous field of vortices may confine spin tex- 
tures in a nearby magnetic semiconductor - thus manip- 
ulation of vortex states can be useful in spintronics [3|; 
the 'ratchet' dynamics of vortices in asymmetric pinning 
profiles is directly related to biological microdevices that 
separate particles by converting random motion into di- 
rected motion [4, 5]; on the microscopic side, the interior 
of a vortex core is fundamentally different from the re- 
mainder of the superconducting sample, and may present 
a unique guiding channel for applications in future elec- 
tronic and optical devices Q. 

One of the most puzzling questions in the area of vortex 
matter in submicron samples is the distinction between 
two allotropies of a vortex state - a 'giant' vortex, where 
all vortices merge into a single singularity, and a multi- 
vortex, where all vortices can be individually resolved. In 
type-II superconductors, transitions between the latter 
two states are of second-order, following the increasing 
lateral compression by e.g. increasing screening currents 
in increasing magnetic field 7, 8\ , or increasing tempera- 
ture which makes the sample effectively smaller in terms 
of the superconducting length scales. Even in numeri- 
cal calculations, it is very difficult to pinpoint the exact 
value of parameters for the giant-to-multi crossover, as 
the order parameter is severely suppressed between vor- 



tices in close proximity. It is therefore no surprise that 
imaging experiments could not verify the existence of a 
giant vortex beyond reasonable doubt |9J. Several years 
ago, Kanda et al. conveyed a clever transport measure- 
ment, where distinction between giant and multi- vortex 
states was made by symmetry matching between the vor- 
tex configuration and the location of several tunnel junc- 
tions [Io[ • Although not always conclusive [22[ , this is the 
best known method to date for giant-vortex detection. 

In this Letter, we present a universal method for the 
observation of formation and decay of multi-quanta vor- 
tex states. Our theoretical simulations indicate that 
the experimentally measured heat capacity of a meso- 
scopic superconductor as a function of magnetic field or 
temperature can unambiguously reveal such transitions. 
The underlying reason can be traced back to the behav- 
ior of the local density of states for quasiparticles, and 
we demonstrate a direct link between the heat capac- 
ity and the sample magnetization. With recent advances 
in calorimetry [llj and magnetometry 12| of submicron 
samples, our findings are of immediate relevance to cur- 
rent experimental efforts. 

The Ginzburg-Landau (GL) formalism has been exten- 
sively used in the past to gain theoretical insight in the 
physics of mesoscopic superconductors. The core of the 
approach is the GL energy functional 
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describing the difference in Gibbs free energy between 
the superconducting (S) and normal (N) state in units 
of Qq = iJ 2 /87r. Here k denotes the GL parameter 
and determines screening of the applied magnetic field 
H from the given superconducting material. In Eq. ([1]) 
all distances are scaled by the coherence length £, the 
vector potential A by cfi./2e£, the magnetic field h by 
H C 2 = ch/2e^ 2 = k\/2H c , and the order parameter ip 
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by its equilibrium value in the absence of the magnetic 
field. The minimization of Q is numerically equivalent 
to solving two coupled GL equations, and for details 
of this procedure we refer to Ref. [13]. Once a sta- 
ble solution is found, we are able to calculate the spe- 
cific heat of the superconducting state from the relation 
C = -Td 2 G/dT 2 Q, as a difference between the total 
heat capacity and that of the sample in the normal state, 
in units of C = H%(0)V/(8irT c ). We start from the equi- 
librium states, increase/decrease the temperature of the 
system by 10~ 4 T C , and calculate numerically the second 
derivative. In what follows, we apply this method for 
a superconducting disk, a simple geometry already ac- 
cessible both theoretically 0] and experimentally (j- 12|. 



Figure [T] shows the energy of all the vortex states found 
in an Aluminum superconducting disk of radius R = 850 
nm and thickness d = 100 nm, at T = 1.1 K (we use 
£(Q) = 100 nm, k = 1.2, and T c = 1.38 K [15]). The con- 
finement effects are more significant for increased vor- 
ticity in increasing field, due to the interaction of the 
flux quanta with lateral boundaries. For that reason, all 
states in Fig. [T] with L > 5 are giant vortices. How- 
ever, for 2 < L < 5 multivortex states can be found at 
lower magnetic field, which are compressed into a giant- 
vortex at higher applied field. This is a gradual, second- 
order transition, and is therefore invisible in the free en- 
ergy curves 0] ■ For clarity, we made distinction between 
multi- and giant-vortex in Fig. Q] by dashed and solid 
lines respectively. In what follows, we discuss the reper- 
cussions of the latter transition on the heat capacity of 
the sample. 

Using attoJoule calorimetry, Ong et al. [16| studied 
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FIG. 1: (color online) The energy of different vortex states 
for an Al superconducting disk of radius R — 850 nm and 
thickness d = 100 nm, at T = 0.8T C , for taken £(0) = 100 
nm and k = 1.2. Solid lines indicate giant- vortex states and 
dashed lines represent multi-vortex states. Insets show the 
density of the superconducting condensate for a L = 2 vortex 
state in multi and giant form. 



the heat capacity of mesoscopic disks as a function of 
the magnetic field and found that the heat capacity is 
directly linked to the vorticity, exhibiting jumps at tran- 
sitions between vortex states. We argue here that the 
heat capacity depends not only on the number of vortices 
in the sample, but also on their configuration. Namely, 
the susceptibility of the sample to heating is linked to 
the kinetic energy of the Cooper-pairs in and around the 
vortex core(s), and the changes in their trajectory upon 
the multi-to-giant vortex transition. Using the definition 
of heat capacity as a second derivative of the Gibbs free 
energy with temperature, we calculate it as a function of 
applied magnetic field for vortex states with vorticity 2 
and 3, both exhibiting multi-to- giant vortex transition in 
Fig. [T] As shown in Fig. [5] for both cases, the general 
trend of increasing heat capacity with field is interrupted 
exactly at the multi-to-giant vortex transition, where a 
sharp decrease of heat capacity is found. 

In what follows, we show that the cause of the observed 
change in heat capacity during the merging of vortices is 
the changing local density of states (LDOS) for quasipar- 
ticle excitations inside the vortex cores. Actually, already 
from the early theoretical works (see Ref. [lj}), we know 
that the bound state spectrum inside the vortex is also 
a function of momentum along the vortex line; as a re- 
sult, the lowest bound state energy for winding number 
L > 1 is L times larger than that of winding number 
1. Therefore, one expects that the low-energy states are 
pushed toward higher energies during merging of individ- 
ual vortices into a giant vortex. To give a quantitative 
measure of this process, we first obtain the order param- 
eter ip and vector potential A of the equilibrium states 
from the GL calculation, which then serve as inputs for 
the microscopic Eilenberger equation fl8| 
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FIG. 2: (color online) The heat capacity as a function of 
magnetic field, for states with vorticity 2 and 3 of the sample 
considered in Fig. [T] Insets depict the vortex configuration 
before and after the multi-to-giant vortex transition. 
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FIG. 3: (color online) The integrated zero-energy density of 
states (LDOS) as a function of the magnetic field for L = 3, 
for a superconducting disk with same parameters as in Fig. 
[1] (a-c) are the representative contourplots of LDOS in the 
disk, at indicated magnetic fields. Interestingly, N(0)[a] fa 
JV(0)[b]«JV(0)[c]. 
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where functions a and b now 



satisfy the independent nonlinear Ricatti equations [1 
In the next step, the LDOS is evaluated from 
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where ry(> 0) is a small real constant. To obtain g(ie„ — > 
E+irj, r, 0), we solve the Eilenberger equations for rj — iE 
instead of the Matsubara frequency u> n . In order to 
find the LDOS, the above equations should be solved for 
a bundle of trajectories with different angle running 
through the given point r and energy e. In our calcula- 
tion we consider only specular reflection for trajectories 
encountering the outer boundary of the sample. 

Early theoretical works already considered LDOS of a 
single vortex [2(J. Further, in Ref. [6] LDOS was cal- 
culated semiclassically for a multivortex vs. the case of 
a giant vortex, for selected vorticities and assumed size 
and distribution of the vortex cores. Here we present the 
full evolution of the LDOS of quasiparticle excitations 
in a mesoscopic superconducting disk, during the multi- 
to-giant vortex transition as a function of the magnetic 
field, where at each step we calculate the distribution of 
the superconducting order parameter. In Fig. [3J we plot 
the zero-energy density of states N(E = 0,T) integrated 
over the sample as a function of the applied magnetic 
field, for the L = 3 vortex configuration. N(E = 0, T) 
increases with applied magnetic field, as in the case of 
an isolated vortex [2(|. When the giant vortex is assem- 
bled from the multi-vortex molecule, the LDOS profile 
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FIG. 4: (color online) The multi-to-giant vortex transition 
revealed through the sharp change in magnetic susceptibility 
as a function of applied magnetic field, showing direct corre- 
lations with the heat capacity, for states with vorticity 2 and 
3. Shaded areas indicate the observed regions of giant-vortex 
formation in both quantities. 



changes from several individual peaks located at each 
vortex to a ring- like bound state with/ without an en- 
closed peak for odd/even vorticity (as in Ref. @). Bound 
states are also found near the sample boundary, due to 
there lowered (non-zero) gap in the presence of strong 
circular Meissner currents. The representative contour 
plots of LDOS for L = 3 are shown in Fig. [3] as in- 
sets. As our main observation, we point out a clear drop 
of N(E = 0) vs. the magnetic field at the multi-to- 
giant vortex transition (see Fig. [3]), where the LDOS 
profile goes through a change of symmetry from three- 
fold to circular symmetric one. We thus confirm that the 
evolution of LDOS with magnetic field is directly linked 
to the specific heat and the thermal conductivity of the 
sample. To enforce this argument, the LDOS can be 
expressed as N(E,T)/N = N(E = 0,T) + a E \E\/A , 
and N(E = 0,T)/N o = i + a k T/T c . Through the re- 
lation C(T)/T = (2/T) dE[E N(E, T)df(E, T)/dT], 
using the Fermi distribution function f(E,T), one ob- 
tains C{T)/{ ln T) - N{E = 0,T)/N Q + a E \E\/A which 
unambiguously shows the link between the calculated 
curves in Figs. [5] and [3] 

The jump of heat capacity between different vortex 
phases can also be expressed using other thermodynamic 
arguments. The discontinuity in the specific heat at a 
phase transition (at field H*) can be calculated as [2l[ 
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where M denotes sample magnetization. Here we apply 
above expression to the multi-to-giant vortex transition, 
where i represents the vortex state of vorticity L just 
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prior, and j represents the L vortex state just after the 
transition. Knowing the result for heat capacity vs. H 
at the multi-to-giant vortex transition, we therefore ex- 
pect to see similar features in the magnetic susceptibility 
X = dM/dH. We calculate the magnetization M as ex- 
pelled magnetic field from the sample M = ((h) —H) / 47r, 
where (h) is the local magnetic field averaged over the 
sample volume. The results of this calculation are shown 
in Fig. |4] They (i) confirm the link between (indepen- 
dently calculated) sharp changes in heat capacity and x 
as a function of the magnetic field, and (ii) show that 
assembly of a giant vortex in superconductors can be de- 
tected even by conventional magnetometry. 

The above prediction is of immediate relevance to ex- 
periments, since both calorimetry and magnetometry are 
readily performed on mesoscopic superconductors. Of 
course, the question of sensitivity and resolution of the 
measurement is an open one, and we address this issue 
in Fig. [5] First, we determined the multi-to-giant vortex 
transition field H* as a function of the size of the Al disk. 
For all considered vorticities (L = 2 — 4), H* was found to 
increase with the radius of the sample. We then scanned 
the heat capacity and magnetic susceptibility versus ap- 
plied field for every size of the sample, and recorded the 
size of the observed jump between values prior and after 
H* . In Fig. [5] we show the absolute and relative size 
of the jump of both magnetic susceptibility (b) and heat 
capacity (c) at temperature 1.1 K. We found that the 
susceptibility shows a clearer signal at the multi-to-giant 



transition for lower vorticity, whereas corresponding dis- 
continuity of heat capacity is more pronounced at higher 
vorticity. Note however that AC and Ax should be di- 
rectly proportional, according to Eq. (U). They indeed 
are, when susceptibility is calculated by x = d 2 Q/dH 2 1 
while we here used the experimental definition of mag- 
netization as the flux expelled from the sample (which is 
directly measured by Hall magnetometry). 

In summary, we demonstrated that second-order tran- 
sitions between multi- and giant-vortex states in meso- 
scopic superconductors can be detected using calorime- 
try. The local density of states for quasiparticles in and 
around vortex cores changes when the vortex configura- 
tion changes, which affects the heating properties of the 
system. The observed sharp change in the heat capacity 
at the multi-to-giant vortex transition can also be linked 
to the magnetic susceptibility, enabling the observation 
of this transition by Hall magnetometry. Our results are 
therefore of immediate relevance to experimental efforts 
in the field, and further work is needed to generalize our 
findings to other systems, such as e.g. Bose-Einstein con- 
densates [l|. 
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FIG. 5: (color online) The absolute and relative size of the 
jump at the multi-to-giant vortex transition, in (a) x(^) 
curves, and (b) C{H) curves, for vorticity L =2,3, and 4. 
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